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ABSTRACT 

DNS studies of dispersed particle motion in isotropic 
homogeneous turbulence [T] have revealed the existence 
of a component of random uncorrelated motion (RUM) 
dependent on the particle inertia t p (normalised particle 
response time or Stoke number) . This paper reports the 
presence of RUM in a simple linear random smoothly 
varying flow field of counter rotating vortices where the 
two-particle velocity correlation was measured as a func- 
tion of spatial separation. Values of the correlation less 
than one for zero separation indicated the presence of 
RUM. In terms of Stokes number, the motion of the 
particles in one direction corresponds to either a heav- 
ily damped (r p < 0.25) or lightly damped (r p > 0.25) 
harmonic oscillator. In the lightly damped case the 
particles overshoot the stagnation lines of the flow and 
are projected from one vortex to another (the so-called 
sling-shot effect). It is shown that RUM occurs only 
when T p > 0.25, increasing monotonically with increas- 
ing Stokes number. Calculations of the particle pair spa- 
tial distribution function show that equilibrium of the 
particle concentration field is never reached, the con- 
centration at zero separation increasing monotonically 
with time. This is consistent with the calculated neg- 
ative values of the average Liapounov exponent (finite 
compressibility) of the particle velocity field. 



1 INTRODUCTION 

Turbulent structures play a crucial role in many parti- 
cle/fluid processes in the environment and industry: in 
powder production, combustion, and the formation and 
growth of PM10 particulates in the atmosphere. An area 
of much investigation is the mechanism for warm-rain 
initiation and in particular the way droplet interaction 
with the small scales of turbulence in clouds influence 
the droplet size distribution. The current consensus is 
that the turbulence demixes the particles leading to a 



much enhanced collision rate (i.e. much greater than 
that based on passive scalar motion [2], [3] )■ Our only 
particular interest in this subject is motivated by the 
processes of turbulent agglomeration and deagglomer- 
ation of a aerosols released in a steam generator tube 
rupture in a PWR. Early experiments and simulations [1] 
have shown that this demixing reaches a maximum when 
the particle response time is typical of the timescale of 
the turbulent structure (i.e particle Stokes numbers ~ 
1), the suspended particles being observed to segregate 
into regions of high straining rate in between the regions 
of vorticity. In addition Maxey and his co- workers [5],[B] 
showed that the gravitational settling of particles in 
homogeneous turbulence was enhanced due preferential 
sweeping in the direction of gravity as particles inter- 
weave through turbulent structures in the flow. Since 
then there have been numerous studies to understand 
and quantify this segregation process. Of particular note 
have been the seminal studies by Collins et al.[7] and 
Wang et al.[8] to quantify the influence of segregation 
on relative dispersion (two particle dispersion) and on 
particle agglomeration and the recent work of Bee [5] 
who expressed the clustering of particle in terms of its 
fractal dimension and showed how this was related to 
the Liapounov exponents of the particle phase space dis- 
tribution. Even more recently Vassilicos[TD] has shown 
that the clustering is strongly related to the acceleration 
stagnation points in the flow as they are swept by the 
large scale motions of the turbulence. 

This paper is not directly concerned with the clustering 
but with an intrinsic property of the motion of iner- 
tial particles in flow fields that are spatially random but 
smoothly varying. Simonin et al. jT| have observed that 
the spatial velocity field resulting from the motion of 
suspended particles in DNS isotropic homogeneous tur- 
bulence consists of two components: a smoothly (contin- 
uous) velocity field that accounts for all particle-particle 
and fluid-particle two point spatial correlations (they re- 
ferred to this components as the mesoscopic Eulerian 



1 Simonin et al. refer to this component as the quasi-Brownian velocity distribution because of its similarity to the spatially uncorrelated 
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(a) Symmetric straining rate flow 



(b) Pairs of counter-rotating vortices 






a) Generation from symmetric strain rate flow 

b) Carrier flow velocity and rms velocity in x-direction 
Figure 1 : Generation of carrier flow field velocity and its x~ variation within a flow lattice cell 



particle velocity field (MEPVF)); and a spatially uncor- 
rected component which we will refer to here as 
(the component of random uncorrelated motion) whose 
contribution to the particle kinetic energy increases as 
the particle inertia increased. Simonin et al [I] at- 
tribute this feature to the ability of the particles with 
inertia to retain their memory of their interaction with 
very distant, and statistically independent eddies in the 
flow field. The flow field in their DNS of homogeneous 
isotropic turbulence as with real turbulence is compli- 
cated with many scales of spatial and temporal exis- 
tence. In our study, the results of which we present here, 
we show the presence of RUM in a much simpler flow 
field of counter rotating vortices and in so doing are able 
to examine the properties of RUM in relation to the par- 
ticle interactions with simple structures. The flow field 
we have used has been used before in studies of particle 
clustering 11 and in particular with the measurement 



of the particle compressibility and intermittency. Our 
eventual aim is to show how these properties are related 
to the occurrence of RUM 



2 DESCRIPTION OF CARRIER FLOW 
FIELD 

This flow field and its properties we have described in 
detail elsewhere [IT] but for completeness we shall in- 
clude a brief description here. We consider dispersion in 
a simple homogeneous turbulent flow field composed of 
pairs of counter rotating vortices which are periodic in 
both the x,y directions with the same periodicity. Each 
lattice cell (the basic periodic element) contains a pair 
of counter-rotating vortices in both the x,y orthogonal 
directions and is constructed from a linear symmetric 
straining flow field in the manner shown in Figure [JJ 



So starting from an initial symmetric straining flow pat- 
tern of width 2L in both the x, y directions (see Fig. 
[JJa) ) , this pattern is repeated front to back in both the 
x, y directions with a strain rate S drawn from a uni- 
form distribution [0, So]- We note that each quadrant of 
this straining rate pattern in FigJTJa) is a quadrant of 
one of the two pairs of counter-rotating vortices formed 
within the lattice cell in FigOTb). As shown in FigOJa), 
the flow velocity u x in the a;— direction has a linear saw- 
tooth profile U(x), with a slope of constant magnitude 
S but with a change in sign across the y-centre line of 
a vortex where the maximum and minimum values 
±SL/2 of U(x) are located: across the x-centre line, 
u x changes to —U(x) as shown in Fig[JJb), consistent 
with the change in direction of the streamlines shown 



in FigUJa). The flow velocity u y in the y-direction at 
(x,y) is —U(y) to preserve continuity of flow through 
out. This cellular flow pattern of counter- rotating vor- 
tices so formed, persists for a fixed life-time selected from 
an exponential distribution with a decay time of Sq 1 , at 
the end of which time, a fresh flow field is generated 
with new values of the life-time and S and the origin 
of the pattern at the same time shifted by random dis- 
placements in both the x and y-directions, drawn in- 
dependently from a uniform distribution [0,2L]. This 
makes the average flow homogeneous with zero mean in 
the x— ^-directions. The important feature of this ran- 
domized flow field is that the equations of motion of an 
individual particle in both the x, y-directions are linear 
and independent of one another (other than through the 



velocity distribution of Brownian particles 

2 the line running in the y-direction passing through the centre of the vortex 



2 



maximum length of time a particle can experience a par- 
ticular value of the straining of the flow in either the x- 
or y-directions before it changes sign). With respect to 
the centre (stagnation point) of a symmetric straining 
flow pattern (see FigOJa)), the flow velocity within that 
flow region is given by : 

u x = +Sx ; u y = —Sy (-L < x < L, —L < y < L) 

(2.1) 

The flow field so generated turns out to be homogeneous 
and stationary but not isotropic. It has the interesting 
property that the Lagrangian fluid point rms velocity 
(along its trajectory) is different from its value at a fixed 
point (Eulerian). 

2.1 Particle equations of motion 

Based on Stokes drag, the particle equation of motion is 

x l + T~ 1 x l + (~iy +1 T- 1 Sx l = 0(2.2) 
— L < Xi < L (i = 1, 2)(xi = x, X2 = y) 



where Xi is measured from a stagnation point. For con- 
venience we normalize Xi on L and express the particle 
response time t p and strain rate S in units of Sq 1 so r p 
is the particle Stokes number and in this case S is drawn 
from a uniform distribution p [0, 1]. We note that when 
the strain rate in the x, y directions the x, y equations 
are those of a damped simple harmonic oscillator and 
as a consequence there are two types of motion, namely 
heavily damped for t p < 0.25 and lightly damped for 
T p > 0.25. This has some bearing on the motion of 
the particle within these vortices. For heavily damped 
motion, a particle remains trapped within a vortex but 
approaches the extremity (the stagnation region) in a 
manner which decreases exponentially with time. On 
the other hand for lightly damped motion the particle 
can escape the vortex or overshoot into an adjacent vor- 
tex, and in so doing overshoot into the adjacent vortex 
to that. This has been referred to as the sling shot ef- 
fect and it is this type of motion that gives rise to RUM. 
Both these types of behaviour are illustrated in Figure [5] 
(a) and (b) for particle response times r p = 0.1 (heavily 
damped) and t p = 1 (lightly damped). 




Figure 2: Particle trajectories in vortices: r p =0.1 , heavily damped; t p = 1, lightly damped 



Figure(2] also show the corresponding values of compo- 
nents Jn and J22 of the the unit deformation tensor J 
given by 9x(xo,i)/cbco where Xo is the position of the 
particle at some initial time say t = 0. In this case we 
have set J%j(0) — 5ij. The equations of motion for Jij 
are obtained from the particle equations of motion by 
partial differentiation wrt xq^ giving 

4 + r^Jij + (-!)*+ = (2.3) 



for which we choose the initial conditions = 
dui/dxj — (—ly^Tp^SSij. This means that together 
with the initial conditions on Jij that J^- (t) = for i ^= j 
and that 

J(t) = |J| = J n J 22 (2.4) 
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a) St = 0.1, heavilydamped 
Figure 3: Compenents of the deformation tensor J 

Note that in this linear system the equation of motion for 
J(t) is the same as that for xi and that the components 
of J are dependent on the position of the particle only 
through the time at which the strain rate experienced 
by the particle changes sign. Fig[3Ja) and (b) shows the 
corresponding values of J(t) along the heavily damped 
and lightly damped (b) trajectories. In both cases the 
value of J approaches zero as t — > oo. However in the 
lightly damped case J(t) passes through zero at inter- 
mediate times as the particle oscillates backwards and 
forwards across a stagnation line. In so doing the value 
of J oscillates from +ve to —ve, with the correspond- 
ing elemental volume rotating through 180° as it passes 
through zero volume. Each time J(t) passes through 
zero, the corresponding particle concentration becomes 
infinite instantaneously. This raises the possibility that 
such events may occur in real turbulent flows and that 
the process of particle dispersion could be a highly inter- 
mittent process associated with large deviations in the 
particle concentrations. It also indicates the conditions 
under which particle trajectories cross. That is given an 
initially uniform distribution of particle with velocities 
uniquely defines at each point by the local fluid velocity, 
as time elapses the particle velocity field is not uniquely 
defines at each point. This ultimately is the origin of 
RUM. 

3 Calculation of RUM 

We calculated RUM by measuring the contribution it 
makes to the particle kinetic energy (per unit mass) . We 
did this by calculating the particle pair velocity corre- 
lation as a function of particle separation. In particular 
we calculated 

R xx (Ax) = (v x 1 \x,0)v x 2 \x + Ax. / 0))/(vl) (3.5) 
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b) St — 1., lightly damped 
and compression \J\ along particle trajectories : 

where (l)and (2)refer to two particles labelled (1) and 
(2)separated by a distance [Ax,0]and (}is an average 
over all particle pairs for which the separation is Ax. 
We began with a fully mixed flow of particles with the 
local fluid velocity. Periodicity of the lattice cell dimen- 
sions means that initially for any integers n, and m,the 
initial particle velocity field satisfies the periodic bound- 
ary condition, v(x, y) = v(x + An, y + Am). Given this 
initial velocity field means that for any given particle 
trajectory Xi(xi,t) there will be an identical particle 
trajectory ~K.2{x\ + An, y + 4mi), for which 

X 2 (xi + 4n, y + Ami, t) = Xi(xi,i) + T; T = (An, Am) 

(3-6) 

That is two particle starting out from identical positions 
within any two lattice cells will experience the same rel- 
ative displacement with respect to their initial positions 
and necessarily experience the same flow velocity at any 
given time even though the whole lattice is displaced 
randomly at each life time. This means that if we con- 
sider particles entering or leaving a fixed cell of lattice 
dimensions, any particle leaving the cell will be replaced 
by an identical particle at the opposite face whose veloc- 
ity is identical to that of the particle leaving. Thus at 
any instant of time, the number of particles within the 
cell remains constant; only the concentration within the 
cell will change. Periodicity of the particle motion and 
homogeneity means further more that for the particle 
pair velocity correlation 

R xx (Ax) = R xx (-Ax) = R XX (A - Ax) = R yy (Ay) 
= R yv (-Ay) = R(A-Ay) 

Because the particle concentration is initially uniform, 
the particle pair correlation is identical to the Eulerian 
spatial velocity correlation which we can calculate di- 
rectly, namely 
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R xx {Ax) 



3 r 3 



I lattice 

After some labour, one arrives at the result 
< Aa; < 2 



u x (x, y)u x {x + Ax, y)dxdy/ (u x ) = — y J u x (x, y)u x (x + Ax, y)dxdy 



R xx {Ax) = [1 - | Ax 2 + ±Ax 3 ] 



2 < Aa; < 4 R xx {Ax) = ^ [12 - l&Ax + 6Ax 2 - §Aa; 3 ] = [9 - 12Ax + fAa; 2 - ±Aa; 3 ] 



(3.8) 



(3.9) 



These forms provide an independent check on our method of calculating the particle pair velocity correlation based 
on counting the number of particle pairs with a given separation . The method within the statistical error correctly 
reproduces the correlation coefficient for the particle initial uniform distribution (See Fig. [5]) 




a) St=,Q.l, time =10 





b) St=0.1, time=40 
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e) St=10, timc=10 

Figure 4: Particle distributions at time t 



f) St=10, timc=40 
10 and 40, for St = 0.1,1., 10 
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Figure 0] shows the positions of 500 particles initially 
distributed uniformally within a cell whose boundaries 
are coincident with a lattice cell initially but which re- 
mains fixed throughout the simulation, and for which 
periodic boundary conditions are applied at the edges of 
the cell to particles entering and leaving it. The result- 
ing particle distributions in Figured] ar e for a given ran- 
dom sequence of the straining rateS'and lattice life time. 
The particular distributions shown are at times t = 10 
and t = 40 for particle response times (Stokes nos.) of 
t p = 0.1, 1.0, 10. We note first that the case of the in- 
termediate response time t p = 1, exhibits the greatest 
degree of segregation, being markedly greater than the 
other two cases. Secondly the segregation increases with 
time and carries on indefinitely. Thirdly there is no iden- 
tifiable self preserving pattern of segregation, nor does 
the segregation align with the local stagnation lines of 
the flow at any instant of time. In fact another ran- 
dom sequence of strain rates and lifetimes would lead 
to an entirely different pattern. The increasing segre- 
gation would lead to a lower fractal dimension and the 
increase of the particle pair concentration at zero separa- 
tion values for the pair separation distribution function 
as shown in Figure |4] described below. 

In calculating the correlation coefficient R xx (Ax)ior 
particle pair separation in the x-direction, the 
positions of N particles within the domain 
[0 < Ax < 4 , < y < 4] are evaluated as a function of 
time by first solving the set of equations of motion 
(|2~2|) . The domain [0 < Ax < 4 , < y < 4] is divided 
up into 10 2 x 10 2 bins of 0.04 dimension in both the 
Ax, y directions, the (m, n)bin being the bin whose 
centre is at Ax =0.04(m — \),y = 0.04(n — \) where 
m,n = 1,..100. For each particle pair whose y- co- 
ordinates satisfy 0.04(n — 1) < y < n0.04 we evaluate 
the separation Ax. If Ax > we determine the value 
of rn for which 0.04(m — 1) < Ax < m0.04 and store 
the value v x 1 \x 1 0)v x 2 \x + Ax,0) in the (m,n)bin. If 
Ax < 0, we store these values in the (m, n)bin for which 
0.04m < 4 - Ax < (m + 1)0.04. In this way, for all 
particle pairs whose x— position is —2 < x < 2 we cover 
the range < Ax < 4. After sorting all particle pairs 
into suitable bins, we then repeat this procedure for M 
realisations of the carrier flow and then sum the values 
of v^^x , 0)v x 2 \x + Ax,0)for all particle pairs within 
each bin and divide this by the number of particle pairs 

to give the average value ^vi 1} (x,0)vi 2) (x -I- Ax,0)Y In 

addition we calculate (v 2 )by summing v 2 over all N par- 
ticles in the box and over M realisations of the flow field 
and then dividing by MN, finally to obtain R xx (Ax) 
using Eq. (|3.5l) . For accuracy versus computation time 
we chose N = 10 3 particles and M = 10 3 realisations of 
the flow. At the same time as summing over the values 
v x 1 \x, 0)vi 2 ^(x + Ax, 0) in all bins with the same value 
of n, we also summed the number of pairs in each bin 
and divided by 0.04, to evaluate the particle pair distri- 
bution g(Ax)which we normalised to unity by dividing 



the number of pairs within the domain < Ax < 4. i.e. 




g(Ax)dAx = 1 



For each of values of r p = 0.1, 1, 10, Figure [5] shows the 
particle pair velocity correlation coefficient R xx (Ax) at 
time t = 20) as a function of the pair separation. These 
three curves are to be compared with the curve for the 
initial distribution given by Eq. p. 91) . It is apparent that 
the case for t p = 0.1 is almost identical to the initial form 
of the pair separation correlation. For the other cases, 
the curves are noticeably different, with the most inter- 
esting property of all, that the extrapolated values for 
Ax = Oare not unity as is the case of the initial pair cor- 
relation function and that of the case of t — 0.1. This is 
most marked for t p = 10. The difference between the ex- 
trapolated value and unity represents the contribution 
of RUM to the particle kinetic energy. More reveal- 
ing are the calculated values of RUM plotted against t p 
shown in FigJBlfor t = 20. We note that there is no real 
measurable contribution from RUM until t p > 0.25. For 
t p < 0.25 all realisable states correspond to a sequence of 
heavily damped system with t p S< 0.25: above t p = 0.25 
the system behaves as a mixture of heavily damped and 
lightly damped harmonic oscillators. The threshold of 
T p = 0.25 for the occurrence of RUM, shows conclusively 
that RUM is associated with the behaviour of a lightly 
damped system in the manner described earlier. 

Finally Figure [7] show the calculated values of the par- 
ticle pair distribution function g{x). Figure[7](a) shows 
that for the 3 values of t p = 0.1, 1.0, 10 the maximum 
peak concentration Ax = Ooccurs at the intermediate 
values of r p = 1.. Fi^JTJb) shows peak distribution in- 
creasing with time for t p = l,a result that is true for all 
values of t p . We recall that these features are consistent 
with the segregation of particles shown in Figure 0] 

4 Summary and Concluding re- 
marks 

We have shown and measured a contribution to the par- 
ticle kinetic energy from a component of random uncor- 
related motion (RUM) of particles suspended in a sim- 
ple random but smoothly varying flow field composed 
of counter rotating vortices. This feature has been ob- 
served in more complex flow. Depending upon the value 
of the particle response time t p (Stokes number) the mo- 
tion of the particles corresponds to that of a heavily or 
lightly damped harmonic oscillator. Within statistical 
error, RUM was only detected when the particle mo- 
tion at some stage corresponded to a lightly damped 
oscillator in particular when t p > 0.25. In this case 
the particles overshoot the stagnation lines in the flow 
and can pass form one vortex to another in the flow. 
This behaviour is consistent with the description given 
by Simonin et al.pQ who first observed RUM in parti- 
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cle motion in DNS isotropic homogeneous turbulence. 
There remains the possibility of linking RUM with the 
occurrence of singularities in the particle concentration 
flow field and to the intermittency in the particle flow 
filed. Simonin et al. [T] have developed a statistical 
approach to describing particle motion in these circum- 



stances. They have considered RUM as a device for 
stabilising the segregation (the RUM acting as a pres- 
sure, the gradient of which acts in opposition to the force 
driving the segregation. Here we note that the segrega- 
tion does not stabilise even though RUM is definitely 
present. 




Figure 5: correlation coefficient R XX (X) at time t = 40, where Xis the ^-separation of particle pairs; note that , 
R xx (X)is symmetric X = 2 
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